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We present a comprehensive list of ground state energies and spin gaps of finite kagome clusters with up to 42 
spins obtained using large-scale exact diagonalization techniques. This represents the current limit of this exact 
approach. For a fixed number of spins N we study several cluster shapes under periodic boundary conditions in 
both directions resulting in a toroidal geometry. The clusters are characterized by their side length and diagonal 
as well as the shortest "Manhattan" diameter of the torii. A finite-size scaling analysis of the ground state 
energy as well as the spin gap is then performed in terms of the shortest toroidal diameter as well as the shortest 
"Manhattan" diameter. The structure of the spin-spin correlations further supports the importance of short loops 
wrapping around the torii. 

PACS numbers: 75.10Jm, 75.10.Kt, 75.40.Mg 



The quest for magnetic materials and model systems ex- 
hibiting quantum spin liquid behavior is of considerable cur- 
rent experimental and theoretical interest. Among such sys- 
tems, the kagome 5 = 1/2 antiferromagnet stands out as a 
prototypical highly frustrated quantum magnet in two spatial 
dimensions 1 that potentially could exhibit spin liquid behav- 
ior. However, a complete understanding of this deceptively 
simple model given by, 

ff = jJ^Si-Sj, (1) 

has proven surprisingly difficult and exact diagonalization 
(ED) results on this model, despite their limitation to very 
modest system sizes, have become crucial for both theoreti- 
cal developments and non-exact numerical techniques. 

Several Exact Diagonalization (ED) simulations were car- 
ried out to explore puzzling facets of the kagome antiferro- 
magnet 214 . ED established some of its remarkable proper- 
ties, such as the absence of magnetic order and the enor- 
mously high number of singlet excitations below the low- 
est spinful excitation. Complementary computational tech- 
niques have also been applied to the kagome Heisenberg 
antiferromagnet including series expansions 15-17 , Quantum 
Monte Carlo 18 , Diagonalizations in the nearest-neighbor va- 
lence bond basis and variants thereof 19-21 , Contractor Renor- 
malization (CORE) 2223 , Multi-scale entanglement renormal- 
ization ansatz (MERA) 24 , and the Density Matrix Renormal- 
ization group (DMRG) 25 ' 26 . 

In this work we report on large scale ED results for the 
ground state energy and the spin gap of the Heisenberg S = 
1/2 antiferromagnet on various kagome samples consisting 
of up to N — 42 spins. This considerable increase in sys- 
tem sizes (and Hilbert space size) was made possible by a dis- 
tributed memory parallelization of our ED codes. We provide 
a finite size scaling of the ground state energy and the spin gap 
as a function of the shortest diameter of the torii, which seems 
to capture the finite size dependence in a more systematic way 
than a simple 1/N scaling used previously 7 ' 8 ' 14 . 
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FIG. 1. (Color online) The two 42-site kagome clusters used in this 
study. The shortest Manhattan diameters, du, are shown as thick 
dashed lines. 



I. THE CLUSTERS 

The list of the clusters used in this study along with their 
properties is given in Table I. For a given number of spins, N, 
we list in most cases several clusters. For a given N there are 
many such clusters and the ones we list are chosen to be close 
to optimal, where optimal refers to a cluster that would have 
the length of both sides as well as the shortest diagonal all 
equal. This is considered optimal since such a cluster would 
be more likely to have the full point group symmetry of the in- 
finite kagome lattice. Of the clusters we consider only 12, 27b 
and 36d have the full symmetry of the kagome plane and they 
are therefore shown in bold. The basis vectors a, b of all clus- 
ters are given in terms of the vectors ai , &2 shown in Fig. 1 
where the two largest clusters 42a and 42b are shown. We 
also list the shortest Manhattan diameter, du, of each cluster. 
This measure can be visualized by tiling the plane with the 
cluster and finding the shortest path between two equivalent 
sites walking along the bonds of the lattice, c^m is then sim- 
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2 


4 


48 


-5.444 875 216 
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-6.589 143 829 
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0.418 800 403 


18 a 
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3 


\/T2 
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-0.448 027 


0.270 115 263 
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48 
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0.284 567 177 


21 
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-9.172 279 619 
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24 
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V7 


\/l2 


6 


96 


-10.589 965 547 


-0.441 249 


0.207 828 742 


27 a 


(2,1), (-3,3) 


V7 


3 


\/l3 


6 


18 


-11.793 996 213 


-0.436 815 


0.275 413 255 


27 b 


(3,0), (0,3) 
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3 


3 


6 


216 


-11.779 504 985 


-0.436 278 


0.268 776 803 


30 
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-0.438 477 
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24 
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-0.438 552 
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48 
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0.170 275 671 
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-15.781 555 118 


-0.438 377 


0.164 189 901 
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sft 
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\/21 
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-17.038 187 797 


-0.436 877 


0.199 163 545 
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\/l3 
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78 


-17.020 192 866 


-0.436 415 


0.222 433 924 
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(-1,3), (5,-1) 
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\/21 


\/28 
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-0.437 999 


0.120 425 


42 b 
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TABLE I. Cluster studied in this work. Listed are: The number of spins N. The basis vectors a, b in terms of ai and a? (each of length 2a). 
The length of the basis vectors, (|a|, b|), in units of 2a along with the length of the diagonal d — min(|a — b|, |a + b|). The Manhattan length, 
d_M, the length of the shortest loop wrapping around the torus. The number of elements of the symmetry group \G\. The total ground-state 
energy, E. The energy per site, E/N. The value of the spin gap, A, between the 5 = ground-state and the lowest S = 1 state for even 
samples or the gap from the S = 1/2 ground-state to the lowest S = 3/2 state for odd samples. Results shown in bold are for clusters with 
the full symmetry of the Kagome plane. 



ply equal to the number of bonds traversed. We also define 
the shortest geometrical diameter, this is simply Min(|a|, |6|). 
We also list the number of elements of the symmetry group, 
\G\, for nearest-neighbor interactions on the cluster. Note that 
clusters 18b, 21, 24, 27b have larger symmetry groups than 
expected based on the applicable symmetries of the infinite 
kagome lattice. 



II. ED ON DISTRIBUTED SYSTEMS 

Typically an ED study of a cluster would involve a complete 
symmetry analysis of the cluster and resulting spectra. One 
could then write the Hamiltonian matrix restricted to a given 
symmetry sector to a file that would be read into memory 
for diagonalization. Using shared memory systems this has 
been achieved in the past for a specific quantum number sec- 
tor for 42 spins on the star lattice 27 . However, this approach 
becomes impossible for the larger system sizes and symmetry 
sectors included in this study, because writing or storing the 
Hamiltonian becomes prohibitively slow and matrix elements 
have to be calculated on-the-fly. Furthermore, the largest sec- 
tors do not fit into accessible shared memory machines any- 
more, such that these large system sizes are only treatable on 
distributed memory systems where each computational node 
only has a relatively small addressable memory space, typi- 
cally a few giga-bytes (or less). In general, it is much easier 



and cheaper to scale a distributed system to a large combined 
memory space than a shared memory system and future large 
ED studies will likely have to be performed on distributed sys- 
tems. Due to the physical constraints of a distributed memory 
system some peculiarities remain: A Lanczos diagonalization 
proceeds by iteratively performing, Na, matrix vector multi- 
plications on vectors of length, M. If treating one element of 
the vector requires a time t, the full cpu-time for the calcula- 
tion is roughly, NuMt. Due to the memory constraints of dis- 
tributed systems the application of a non-trivial symmetry that 
would reduce M by a factor of K almost always will increase 
t by a factor larger than K, even when the number of cpu's 
used is the same. The result is that it is slower to diagonalize 
the smaller symmetry reduced Hamiltonian than the full one. 
For the largest clusters we therefore only use very few simple 
symmetries. For the 42 site clusters when no lattice symme- 
tries were used the resulting maximal Hilbert space dimension 
is then Q /2 = 269'128'937'220. 



III. RESULTS 

Our main result is the table 1 listing the ground state energy 
and the spin gap of the samples considered in this work. Some 
of the smaller samples have already been studied in the past 
and energies were quoted for 12,15,18b, and 21 in Ref. 3 28 , 
for 27b in Refs. 5 and 7 and for 36d in Refs. 5 and 8. An 
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FIG. 2. (Color online) Ground state energy per site plotted as a 
function of the inverse shortest geometrical diameter. Remarkably 
good agreement with recent DMRG results -6 on long cylinders with 
the same diameter is observed. For comparison we also display the 
MERA upper bound 24 , series expansion results for the Marston-Zeng 
valence bond crystal 1617 , as well as the DMRG upper bound 26 . 
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FIG. 3. (Color online) Spin gap of even and odd kagome samples 
obtained by ED and plotted as a function of the inverse shortest geo- 
metrical diameter. For comparison recent DMRG data -6 obtained on 
long cylinders with the same diameter are shown. 



approximate ground state energy and spin gap of the sample 
36c have been obtained by DMRG in Ref. 25. 

We plot the ground state energies as a function of the in- 
verse geometrical diameter, which in our convention corre- 
sponds to l/|a| and display the result in Fig. 2. This presen- 
tation seems to capture the finite size effects in a more sys- 
tematic way than the previously used 1/N scaling, as the data 
seems to behave consistently upon increasing the system size, 
while keeping the diameter constant. Furthermore we observe 
good agreement with recent DMRG data 26 on long cylinders 
with the same diameter, thus corroborating the accuracy of the 
DMRG simulations. 

Next we plot the spin gap data in Fig. 3 in the same way. For 
each diameter we observe that the spin gap is monotonously 
decreasing with system size (for even and odd samples sep- 
arately). While it is difficult to extrapolate the ED spin gap 
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FIG. 4. (Color online) Magnified spin-spin correlations in the ground 
state of the 30 and 36d kagome samples. The 30 site cluster with its 
low symmetry requires three distinct reference sites to be consid- 
ered. The diameter of the circle is proportional to the magnitude of 
the spin-spin correlation with the reference site (black filled square). 
Blue (red) color denotes positive (negative) correlations. The nearest 
neighbor correlations (all antiferromagnetic) have been dropped for 
clarity. One observes pronounced staggered spin-spin correlations 
along selected loops wrapping around the torus in the left two panels 
and the lower right panel. 



for constant diameter due to a lack of system sizes for large 
diameters, the qualitative agreement with DMRG results on 
long cylinders suggests that it is indeed the diameter which 
controls the finite size effects upon moving towards the two- 
dimensional bulk limit. 

Complementary evidence for the important role played by 
short loops wrapping around the torus is provided by a mag- 
nified view on the spin-spin correlations in the ground state of 
several even samples displayed in Fig. 4 (see Ref. 5 for tabu- 
lated values of sample 36d). In a spin liquid with a very short 
correlation length one would expect spin correlations between 
distant sites to be very weak and also not to depend signifi- 
cantly on the sample geometry. Indeed in Fig. 4, most spin- 
spin correlations are quite weak, but pronounced staggered 
spin spin correlations along selected loops wrapping around 
the sample are revealed. We expect these resonances to disap- 
pear once the samples are sufficiently wide. 



IV. CONCLUSIONS 

We have presented a systematic study of the ground-state 
energies and spin gaps of many Kagome clusters. In particu- 
lar we have presented results for N = 39, 42 obtained on dis- 
tributed memory parallel clusters. Future quantum numbered 
resolved ED studies on the largest samples should allow to 
check the appearance of the topological degeneracy required 
for the Z2 spin liquid advocated in Ref 26. 
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